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Abstract 

According to the Svetitsky-Yaffe conjecture, a (d + l)-dimensional 
pure gauge theory undergoing a continuous deconhnement transition is 
in the same universality class as a d-dimensional statistical model with 
order parameter taking values in the center of the gauge group. We show 
that the plaquette operator of the gauge theory is mapped into the energy 
operator of the statistical model. For d = 2, this identification allows us to 
use conformal field theory techniques to evaluate exactly the correlation 
functions of the plaquette operator at the critical point. In particular, we 
can evaluate exactly the plaquette expectation value in presence of static 
sources, which gives some new insight in the structure of the color flux 
tube in mesons and baryons. 
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1 Introduction 


Consider a (d + l)-dimensional pure gauge theory undergoing a continuous 
deconfinement transition at the critical temperature T^,. The effective model 
describing the behavior of Polyakov lines at finite temperature T will be a d- 
dimensional statistical model with global symmetry group coinciding with the 
center of the gauge group. Svetitsky and Yaffe 0 were able to show that this 
effective model has only short-range interactions. If also the d-dimensional 
effective model displays a continuous phase transition then it follows from uni¬ 
versality arguments that it belongs to the same universality class of the original 
gauge model. 

Therefore all the universal properties of the deconfinement transition can be 
predicted to coincide with the ones of the dimensionally reduced effective model. 
These include the values of the critical indices, the finite-size scaling behavior, 
and the correlation functions at criticality. The conjecture has passed several 
numerical tests, which became more and more stringent in the last years due 
to the increased precision reachable with Monte Carlo simulations (see 0 and 
references therein). 

It is clear that the Svetitsky-Yaffe conjecture becomes very predictive for 
d = 2, where, using the methods of conformal field theory, the critical behavior 
can be determined exactly. For example, the critical properties of (2 -|- \)D 
SU{2) gauge theory at the deconfinement temperature coincide with those of 
the two-dimensional Ising model. This allows us not only to predict the exact 
values of the critical indices, but also to write down all the multipoint correlation 
functions of the Polyakov loop at criticality ||] . 

What is needed to fully exploit the predictive power of the Svetitsky-Yaffe 
conjecture is a mapping relating the physical observables of the gauge theory 
to the operators of the dimensionally reduced model. In d = 2 the knowledge 
of this mapping is equivalent to solving the gauge theory at the deconfinement 
temperature. 

The correspondence between Polyakov line and order parameter of the ef¬ 
fective model is the first entry in this mapping and is intrinsically contained 
in the Svetitsky-Yaffe conjecture. It is natural to ask what operator in the d- 
dimensional model corresponds to the plaquette operator of the gauge theory: 
symmetry considerations suggest the energy operator as a natural candidate. 
In this paper we show that this is actually the case, and we describe some 
consequences of this identification. 

The correctness of the identification plaquette-energy is shown in Sec. 2 by 
studying the finite-size behavior of the plaquette operator in (2-|- 1)Z? Z 2 gauge 
theory at the critical temperature. We show that it coincides with the (highly 
non-trivial) finite-size behavior of the energy operator in the 2D Ising model at 
criticality. 

In Sec. 3 we compute correlation functions of the plaquette operator by 
using conformal field theory techniques. In particular, the expectation value of 


1 


the plaquette in the vacuum modified by static sources can be computed for 
(2 + 1)1? SU{2) and SU{3) gauge theories at the deconfinement temperature, 
providing physical insight about the structure of the color flux-tube in mesons 
and baryons. 


2 Finite—size behavior of the plaquette expecta¬ 
tion value 


Finite-size effects at criticality are typically rather strong, due to scale invari¬ 
ance, and non-trivial. Therefore they are ideally suited to compare theoretical 
predictions with, for example, results of Monte Carlo simulations. In particular, 
for two-dimensional statistical systems, the critical behavior, including finite- 
size effects, is completely understood with the methods of conformal field theory 
(CFT). We want to exploit this fact to establish the correspondence between the 
plaquette operator in a (d -I- l)-dimensional lattice gauge theory at the decon¬ 
finement transition and the energy operator of the corresponding d-dimensional 
statistical model. 

Consider for example the 2D Ising model: the shape and size dependence at 
criticality of the expectation value of the internal energy on a torus is given by 

ill 

W = (1, 

where A and r are respectively the area and the modular parameter of the torus, 
and Zi /2 is the Ising partition function at the critical point: 



!/=2 


0AQ,r) 
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( 2 ) 


where are the Jacobi theta functions and rj is the Dedekind function (for 
notations and conventions see Ref.|^). 

Comparing Eq. (i with the finite-size behavior of the plaquette operator in 
a 3D lattice gauge theory such that the center of the gauge group is Z 2 provides 
a stringent test of our identification. The simplest choice is the 3D Z 2 gauge 
model, for which it is possible to achieve very high precision in the Monte Carlo 
evaluation of physical quantities, and accurate estimates of the deconfinement 
temperature are available. 

Therefore we considered the (2 -|- 1)1? Z 2 gauge model on lattices of size 
Li X L2 X Lt where Lt <C ^1,^2 with periodic boundary conditions on all di¬ 
rections and we studied it at the critical coupling /3c (Tt), which is known to 
high accuracy for several values of Lt [Ij. By performing Monte Carlo simula¬ 
tions at different values of Li,L 2 , we can compare the finite-size behavior of 
the plaquette expectation value with Eq. {§■ 
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More precisely, we will show that the plaquette operator is a mixture of the 
identity and energy operators of the 2D CFT: on one hand, both these operators 
transform as singlets under Z 2 , and therefore can contribute to the plaquette 
operator; on the other hand, we know that the plaquette expectation value does 
not vanish in infinite volume, unlike the energy operator of the 2D Ising model 
(see Eq.(|Il)). Therefore, a non-vanishing contribution of the identity operator 
must be expected in the plaquette expectation value. Hence our conjecture is 

(□) =ci(l) + c,(e) (3) 

where the expectation value in the l.h.s. is taken in the LGT, while the ones 
in the r.h.s. refer to the CFT. The prediction for the finite-size behavior of the 
plaquette expectation value is therefore 

(FI)liL2 = Cl + Ce +C>(1 /LiL2) (4) 

where F is a function of the modular parameter only: 

Smr |?7 (t)|^ 

Zu2{r) ^ ^ 

We performed our Monte Carlo simulations at Lj = 6 with critical coupling 

= 6) = 0.746035 (6) 

We measured the plaquette expectation value for lattices of area 100 < 
A < 6400 and asymmetry ratio ‘Amr = 1,2,4. The space-like and time-like 
plaquettes have different expectation values, therefore must be fitted separately 
with Eq.(||). The Monte Carlo simulations were actually performed in the ZD 
Ising (spin) model, which is exactly equivalent through duality to the Z 2 gauge 
model. This choice allowed us to use a non-local cluster simulation algorithm. 

The agreement is very good, giving = 0.7 for space-like plaquettes and 
Xred = 0-9 for time-like plaquettes. These data are plotted in Eig. 1. 


F(r) = 


rV 


3 Correlation functions of the plaquette opera¬ 
tor 

In this section we will exploit the new entry we added to the Svetitsky-Yaffe 
mapping to compute correlation functions of the plaquette operator at the de¬ 
confinement temperature. This will provide some new insight into the structure 
of color flux tubes in mesons and baryons. 

^ It must be noted that the expectation values of space-like and time—like plaquettes for a 
given lattice are not statistically uncorrelated, since they were extracted form the same sample 
of configurations. 
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Figure 1: Size and shape dependence of the plaquette expectation value. Black 
dots correspond to square lattices = \). White dots and squares corre¬ 

spond to rectangular lattices with Qmr = 2 and Qttit = 4 respectively. Both 
time-like and space-like plaquettes are shown, the latter having lower expecta¬ 
tion values. The lines correspond to the best fit to iJg. 
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Consider for example (2 + 1)D SU{2) LGT at the deconfinement temper¬ 
ature. To study the flux tube structure in a “static meson” we can consider 
the plaquette expectation value in the vacuum modified by the presence of two 
static sources, i. e. the correlation function of the plaquette operator with two 
Polyakov loops: 


G{x,Xi,X2) = {a{x)P{xi)P{x 2 )) - {a){P{xi)P{x 2 )) (7) 

where x, xi, X 2 are points in the 2D space. This will be given by the correlation 
of the energy operator with 2 spin operators in the 2D critical Ising model: 


G(a:,a:i,a: 2 ) oc {e{x)a{xi)a{x 2 ))ising (8) 

The r.h.s. is easily computed in CFT and we find 


G(x, Xi,X 2 ) oc - 

(|a; 

We have plotted this function in Fig. 2. 


\X1 - X2f^^ 

-Xi\\x- X2\f^'^ 


( 9 ) 



Figure 2: Structure of the flux tube in a ’’static meson” at the deconfinement 
temperature, Eq. §)■ 
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More interesting is the case of (2 + 1)1? SU (3) LGT, where we can consider a 
’static baryon” by modifying the vacuum with three static sources and compute 

G{x,Xi,X2,X3) = {a{x)P{xi)P{x2)PiX3)) - {O) {P{xi)P{x2)P{x3)) (10) 


Our identification gives 

G(x,Xi,X2,X3) oc (e(a:)cr(a;i)cr(x2)cr(a;3))3 

— state Potts ( 11 ) 


where the correlation function on the r.h.s. must be computed in the c = 4/5 
CFT describing the three-state Potts model at criticality. This is done using 
the methods introduced in Ref.0 (see also [||) and gives 


G{x,Xi,X2,X3) cx 


(kl - X2\ \xi - Xsl \X2 - 
{\x - a;i| Icc - X 2 I \x - X3\f^^^ 




l/i(2/)p 


9r3(3/5)r(i/5) 

4r3(2/5)r(4/5) 


\h{y)Y 


( 12 ) 


where, introducing a complex coordinate z in 2D space, y is the conformally 
invariant cross-ratio 


_ (Z - Zi)(z2 - Z3) 
(Z - Z3)(Z2 - Zi) 


(13) 


and /i and /2 are hypergeometric functions: 


h{y) = +(4/5,7/5; 8/5; y) (14) 

f2{y) = y-3/®+(l/5,4/5;2/5;y) (15) 

G{x,Xi,X 2 ,X 3 ) is plotted in Fig. 3, for the case in which the three static 
sources form an equilateral triangle. Notice that this calculation brings strong 
support to the “Y” structure of flux tube in baryons (see e.g. Q and references 
therein), as opposed to the “A” structure |^. 


4 Conclusions 

In this paper we have added a new entry to the Svetitsky-Yaffe mapping between 
(d+l)-dimensional LGT’s at the deconfinement temperature and d-dimensional 
critical statistical models, namely we have shown that the plaquette operator of 
the LGT is mapped into the energy operator of the statistical model. 

For d = 2, this identification allows in principle the exact evaluation of all 
correlations of the plaquette operator at the deconfinement point, providing a 
useful tool for the study of the color flux tube in mesons and baryons. 
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Figure 3: The structure of the flux tube in a ’’static baryon”, Eq. 
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